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The plane-wave dynamics of Soj — > 2a;, subharmonic optical parametric oscillators containing 
1^-^ [ a second harmonic generator of the idler-wave u> is analyzed analytically by using the meanfield 

approximation and numerically by taking into account the field propagation inside the media. The 
resonant x'^'(~3'^i 2(j-',<^) : x*'^'' (~2a;; tJ, cascaded second-order nonlinearities induce a mutual 
injection- locking of the signal and idler waves that leads to coherent self phase-locking of the pump 
■ and subharmonic waves, freezing the phase diffusion noise. In case of signal-and-idler resonant 

devices, largely detuned sub-threshold states occur due to a subcritical bifurcation, broadening out 
the self-locking frequency range to a few cavity linewidths. 
O ■ PACS: 42.65.-k, 42.65.Yj, 42.65.Sf, 42.65.Ky 
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I. INTRODUCTION 

^ , The synthesis of phase-coherent (sub)harmonic optical radiation {uj,2uj, ...,iVw) is useful in high precision optical 
^ measurements, such as optical frequency metrology [Q, In frequency metrology, a phase-coherent optical by-N 
CO divider allows to reduce the absolute measurement of an optical frequency Nu (e.g, couting cycles in the hundreds of 
THz range) to the measurement of a smaller frequency interval A = w by use of a femtosecond laser comb generator 
whose radio- frequency intermode spacing is phase- locked to a primary microwave clock . The generation of a comb 
of phase- locked harmonic radiation can also be the starting basis for the generation of attosecond pulse train by Fourier 
' synthesis [Q, Indeed, the superposition (Fourier synthesis 0, ||^) of a comb of iV equal amplitude optical 

■~~>~ : harmonic fields with a controllable relative phase relationship can lead to a temporal train of ultra-short pulses with 
^ a sub-optical cycle duration (t ~ 2tt/Nlu) and a repetition rate equal to the fundamental harmonic uj. Phase-locked 
• ^ optical harmonic generators are often based on many lasers linked via nonlinear upconversion parametric processes 
' lH ' [0 ' (ill ' 0] • Phase coherence among the harmonic waves is usually achieved by use of complex electronic 
(~| phase-locking loops. The high conversion efficiencies of optical parametric oscillators (OPOs) combined with the 
f*) i- engineering flexibility offered by periodically poled (PP) nonlinear materials may allow a compact implementation of 
such optical subharmonic generators or attosecond pulse generator. 

The study of OPOs as by-A^ dividers (A^ = 2,3,4) of a pump frequency I'p is motivated by their capacity to 
perform the phase-coherent division of a pump photon into two highly phase-correlated subharmonic photons. In 
H , precision measurement setups, subharmonic generation leads to a subsequent phase noise reduction (-^A^^) compared 
to harmonic generation {xN"^). Graham and Haken have first demonstrated that the phase sum tpi + ip2 of the 
idler and signal waves from an OPO follows adiabatically the phase noise ipp of the pump laser, while their phase 
difference ip2 — fi undergoes a phase diffusion process stemming from the interaction of both modes with the vacuum 
fluctuations [p^ . To act as pertinent phase-coherent dividers, the idler and signal waves have thus to be phase-locked, 
for instance by use of an electronic servo control that forces ip2 — V?i to copy the phase of an external RF oscillator 
referenced to a primary clock. Such an electronic servo technique has been implemented to develop phase-coherent 
divide-by-2 (2:1) |T^, and divide-by-3 (3:1) OPOs for high resolution optical metrology ||l^, with corresponding 
residual phase difference variance well below 1 rad^. 

Electronic phase-locking loops are however subject to limited response time bandwidth and require fast OPO cavity 
length actuators. As an alternative, all-optical self phase-locking (SPL) methods are currently being investigated p7| . 
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| [18[ to by-pass the bandwidth hmitation of electronic servos and simphfy the implementation of phase-locked dividers. 
These methods are based on the self injection-locking of the subharmonic waves, similarly to the injection locking 
process of a slave laser oscillator by a master laser which possesses a better spectral purity and frequency stability , 
||19| . In injection-locked systems, the phase of the slave oscillator, <j){t) = mt + ip{t), is a perfect copy of the phase of 
the injecting master oscillator. Phase locking can occur only when the two frequencies are close enough, e.g. within a 
certain locking range Aujiock — — whose extent is proportional to the squared ratio -y/ Pm/Ps of the injecting 
master power to the slave power. It is thus expected that the self-locking range in an OPO divider should depend 
on whether the OPO is configured as a doubly-resonant (DRO), a triply-resonant (TRO) or a pump-enhanced singly 
resonant oscillator (PRSRO). Strong self injection- locking regime occurs only when the signal and idler waves are 
simultaneously resonant. 

For divide-by-2 OPOs based on a type-I nonlinear process (the signal and idler are identically polarized), self- 
injection locking occurs naturally when the OPO is tuned close enough to the frequency-degeneracy |^^. Recently a 
type-II self phase-locked 2:1 DRO was demonstrated, in which a small mixing of the two orthogonally polarized signal 
and idler waves, performed by an intracavity wave plate, induces the injection locking ||l^. In both cases the evidence 
of self phase locking was provided by the high level of phase-coherence between the frequency degenerate signal and 
idler waves. Following the experiment in ref the theory of linearly coupled type-II 2:1 SPL-DROs has been 
reported |21|. In the case of 3:1 OPO dividers, obviously such a linear coupling cannot be implemented, and one must 
use a nonlinear coupling via x*-^-* (—3a;; 2uj, uj) : x^^-* {—2uj; uj, ui) cascading processes (for instance by introducing inside 
the OPO cavity a second nonlinear material phase-matched for the SHG of the idler wave). The first implementation 
of such a self phase-locked (SPL) OPO was recently reported, and used a dual-grating periodically poled lithium 
niobate chip carrying the OPO and SHG sections in a pump-resonant idler-resonant configuration (PRSRO) |l^. An 
extremely small locking range of ~500 kHz, corresponding to a fraction of the idler cavity linewidth, was reported 
due to the weak (non resonant) doubled idler power. A different situation should occur when both the signal and 
idler waves are resonant, as in a DRO or a TRO. Because of the enhanced energy flow exchange between the 2lu and 
UJ modes, the dynamics of such a self-phase locked OPO (SPL-OPO) is expected to differ significantly from that of a 
conventional OPO. Strictly speaking the 3:1 SPL-DRO or TRO can be regarded as a degenerate sub/second harmonic 
generator (2a; ^ a; -fa;) internally pumped by the signal wave of a non-degenerate OPO (3a; ^ 2a;-|-a;). While theories 
of OPO devices containing additional up or down conversion of the signal and idler waves have been investigated in 
the past with the aim of generatin g ne w frequencies [ p2[ , , , [ p5[ , or revealing the quantum noise signature of 
competing x'^^ nonlinearities |26|, ||27| such a subharmonically resonant configuration of competing nonlinearities in 
OPOs has never been theoretically investigated in detail. A 3:1 DRO-SHG rough rate equation analysis based on 
photon flux conservation, but neglecting the role of cavity detunings and field phase coherence, was previously given 
by Zhang et al with the conclusion that self-phase locking should manifest through an imbalance of the signal and 
idler intensities at exact 3:1 division. Actually our in-depth analysis shows that such an imbalance holds only for zero 
detunings: When the nonlinear phase shift due to cascading is compensated by a non-zero optimal cavity detuning, 
one can obtain the same output intensities as in a conventional DRO. 

In this paper, we provide a detailed plane-wave description of the dynamics of divide-by-3 SPL-OPOs that encom- 
passes all the three cavity configurations (only the main results are summarized in the conclusion for the PRSRO, 
which will be detailed elsewhere). The dynamics of the nonlinearly coupled OPO is shown to differ substancially from 
the linear coupling case of a type-II 2:1 SPL-DRO In particular, the nonlinear coupling gives rise to a subcrital 
bifurcation for any non-zero cavity detuning while conventional OPOs undergo subcritical bifurcations only for the 
case of largely pump-detuned TROs Q , ||3^ . As a consequence of the nonlinear phase shifts due to the cascading 
processes, SPL-DROs or TROs will be shown to display large-intensity stable stationary states that would correspond 
to non-lasing detuning domains in conventional OPOs, extending thus the phase-locking detuning range up to a few 
cavity linewidths. 

This paper is organized as follows. In section II we describe the basic plane-wave ring cavity model leading to the 
meanfield solutions that will be compared with the full propagation solutions. In section III, we first treat in detail 
the DRO configuration case. We then extend the theory to the case of a TRO in section IV. Main results on the 
PRSRO configurations will be also summarized in the conclusion. Section V discusses the practical implementation 
of such SPL-OPOs on the basis of the theoretical findings and section VI concludes with some prospective studies 
aroused by this first theoretical approach. 
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II. BASIC MODEL AND EXACT NUMERICAL SOLUTIONS 



The OPO-SHG device described in the paper is schematically sketched in Fig.l. We consider a ring cavity containing 

the cascaded OPO and SHG nonlinear media, each of length Li and L2 respectively and assume that the pump, signal 
and idler waves satisfy already the 3:2:1 frequency ratios. Though the drawing depicts the case of a dual-grating PP 
material, phase- matched for an (eee) interaction, the model is valid also for any separate birefringent material sections, 
or a single material phase-matched for both interactions. In the case of birefringent phase-matching however, the 
OPO should be of type-II kind (e.g. e ^ o + e oy o ^ e + o) and the SHG of type-I kind (resp. e + e ^ o oy + ^ e), 
so that the system is described by only 3 field variables. 



A. Plane-wave propagation equations 

We denote by 2; = Z/Li and z' = Z/L2 the normalized propagation distances within each crystal, such that 
< .2 < 1 between points "0" and "1" and < 2' < 1 between "1" and "2". Let Ej (in m/V) be the slowly varying 
complex field amplitudes where the subscripts j = p, 2, 1 stand respectively for pump, signal and idler 

ej{Z, t) = ^Ej{Z, t) exp [i{ijOjt - kjZ)] + c.c. (1) 

and let Nj be the complex field variable, such that liVjl^ is the number of photons in mode j at plane z inside the 
ring cavity, by 



N, = .b£}LEj (2) 



where V is the average volume occupied by the modes inside the resonator. Then, the reduced field amplitudes Aj 
are introduced by scaling Nj with the small signal gain coefficient giLi of the OPO crystal 

A, = giLiN, (3) 

with 



^ dppo / 2nUJiUJ2L0p ... 

c y eoVnin2np 

In Eq.(4), dopo = ' [x'^^H^^o'; 2a),ti))/2]: 62 • Oi (in m/V) is the effective nonlinear coefficient for the OPO 
interaction and rij are the refractive indices at frequency cuj . The reduced Maxwell equations for these field amplitudes 
Aj (z,t), when propagating through the (3a; 2u),iv) crystal without any diffraction eflfect, are 

^ = iAiA2 (5a) 

^ = iApAl (5b) 



dz 
dAi 

dz 



iApA^ (5c) 



using the standard slowly varying amplidude approximation and the usual change of variables z — > z, t — > t — nz/ c, 
where Ti is the mean linear refractive index. In Eqs. 5 perfect phase- matching has been assumed, /^kopo — kp — 
k2 — ki = 0. The propagation equations in the SHG crystal, allowing for a non-vanishing wavevector mismatch 
^ksHG = k2- 2fci, are 

^.0 (6a) 

^ = iSAl exp(+i2^2') (6b) 
= iSA2Al exp(-i2^2') (6c) 
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with initial conditions Aj{z = 0) = Aj{z = 1). 

The two cascaded crystals are usually phase-matched by adjusting their temperature or angle. But, exact zero 
mismatch for both interactions may be difficult to achieve, so that a phase-mismatch parameter ^ = AksHGL2/2 is 
introduced. The parameter S in eqs. (6) is the ratio of the SHG to OPO small signal gains 



9iLi 

where the small signal gain SHG coefficient (in m^^) is 



(7) 



dsHG / 2flUjfuJ2 ,„s 
52 = \ 2 8 

For a type-I (eee) PP crystal the effective nonlinear coefHcients satisfy dopo ~ dsHG, so that with ni « n2 ~ Up 
and ujp — 3ijj, UJ2 — 2cj, uJi ~ uj^ the relation (7) can be approximated hy S ~ {L2/ Lij/^/i. For L2/L1 = 1/3 for 
instance, one has S ~ 0.2. 

The solutions at the SHG crystal exit can be limited to few terms of the Mac-Laurin field expansion serie at the cell 
entrance, because of the smallness of the parametric gain coefficients. We shall refer to MLl approximation, keeping 
only the lowest-order perturbative terms, i.e. quadratic terms in powers of the field amplitudes. At the exit of the 
SHG crystal, the MLl axpproximation provides the solutions for the signal and idler amplitudes as functions of their 
value at the entrance of the OPO crystal (point "0" in Fig.l) 

Ap{t, Li + L2) ^ Ap{t, 0) + iAi{t, Q)A2{t, 0) (9a) 
A2{t, Li + L2) - A2{t, 0) + iAp{t, Q)A\{t, 0) + ix*Al{t, 0) (9b) 
A, {t, Li + L2) = A,{t, 0) + lApit, 0)A* {t, 0) + ixA2it, 0)AI {t, 0) (9c) 

with the nonlinear coupling constant x 

X = 5exp(-*0(sin^/0 (10) 



B. Boundary conditions 



The reduced Maxwell equations have to be completed by the boundary conditions for the three waves at the entrance 
of the OPO crystal to derive the cavity equations. The field Aj{t + r, 0) at location "0" and at a time t + t , where 
T = A/c is the cavity roundtrip time (A is the total cavity optical path), is the field Aj{t,Li + L2) at point "2" of 
Fig.l, which propagates freely after bounces at the totally reflecting (i? — 1) and output coupling {rj,tj) mirrors and 
eventually summed with an input field, as it is the case for the pump field. We define by rj (tj) the overall mirror 
amplitude reflectivities (transmissivities) , such that r| + 1| = 1 and define the amplitude loss coefficients kj ^ 1 — rj . 

Then, the boundary conditions take then the form 

Ap{t + T,0) = rpexp{iAp)Ap{t, Li + L2) + A,n (11a) 
A2it + T,0)= r2 exp(iA2)A2(i, Li + L2) (lib) 
Ai{t + r, 0) = ri exp(iAi)Ai(t, Li + L2) (11c) 

where Ain denotes the input pump amplitude inside the cavity, Ain — tpAf^^ . Each Aj is the usual cavity detuning 
between the frequency of the waves and the corresponding cold cavity frequency, scaled to the HWHM cavity resonance 
width ( Only nearly resonant waves will be considered, Aj << 2tt). 

The set of equations (9) and (11) forms a mapping of the field amplitudes at location "0", from which meanfield 
equations may be derived. 

The meanfield model was originally derived from a similar set of equations for a two-level cell ring cavity device when 
the atomic dephasing time is much greater than the round-trip time r, but much smaller than the photon lifetime t/k 
. For an OPO, the situation is somewhat different because the response of the crystal is assumed to be instantaneous 
(see Eqs. 5). Nevertheless, meanfield equations can be also derived from the MLl equations (9) and the boun dary 
conditions (11) if the field amplitudes are slowly varying during a round-trip time, i.e. T{dAj/dt) << Aj{t) 
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In the conventional TRO case (|x| =0 , Aj << l,Kj << 1), nontrivial homogeneous stationary solutions exist for 
the idler and the signal only if the following relation holds , |^ 

A= Ai/ki = A2/K2. (12a) 

and the intracavity pump intensity Ip — \Ap\^ is clamped to the threshold input intensity, 

Ith = + A1A2 (12b) 

whatever the input intensity might be above the threshold. In the case of the DRO, the relation (12a) still 

holds. But solutions for the signal and idler intensities can be determined only if the expansion (9) is continued up 
to the cubic terms, leading to 



h = 2k2 



+ A)ilp/lth-l)-l 



(12c) 



J1//2 = A2/A1 = (12d) 
cot{ipp - (fii - 1P2) ^ - A (12e) 

The bifurcation is supercritical for any detuning . Besides, the phase difference tpi — ip2 is an undetermined quantity 
while the sum phase (fi + ip2 depends on the pump laser phase pl| . Actually, this classical phase indetermination is 
compatible with the result of the quantum fluctuation theory of parametric oscillators ||l^ . 

III. SPL-DRO CASE 
A. Stationary solutions 

In the case of a DRO, only the signal and idler waves resonate with the cavity while the pump wave is a travelling 
wave {tp = 1, Tp — 0). Assuming ki^2 <C 1, Ai,2 1, Eqs.(9-ll) provide the stationary solutions, Aj{t + r, 0) = 
Aj{t,0) = Aj, which are also those of the meanfield model, 

Ap — Ain (13a) 
(k2 - iA2)A2 = %ApA\ + ix'Al (13b) 
{Kx-iAx)Ax=iApA2 + ixA2A\ (13c) 

Eqs (13) can be conveniently solved by setting Aj — aj ex.p(iipj), giving rise to 

{k2 — iA2)a2 — iaptti exp(z/i) + i \x\ al exp(— 277) (14a) 

{ki — iAi)ai = iapa2 exp(i/i) + i \x\ 0L\a2 exp(+i77) (14b) 

[i = ipp-Lp2-V\ (14c) 

77 = (/?2 - 2v3i - ^ (14d) 

The non-trivial solutions (^Aj ^ 0) can be easily handled by setting X = (|x| ai)^ and Ip — a^. Then the scaled 
idler intensity is found to be solution of 

- 2{Ip - KiK2 + AiA2)A + |(ki - iAi){K2 + iA2) - Ip\^ = (15) 

Let us note that the coupling terms in Eqs. (14) allow to obtain the non-trivial solutions at the MLl approximation 
order, unlike for the conventional DRO. 

Three cases have to be distinguished in order to solve Eq. (15), Ai^2 — 0, A1A2 ^ (but with A1A2 > 0), and 
Ai,2 = 0, A24 0. 

For A1A2 7^ 0, there are two solutions 



X± = (Ixl ai)\ =Ip- KiK2 + A1A2 ± 2^AiA2 [Ip - lo] (16) 
for any input Ip > Iq, where Iq is the input intensity at the saddle- node bifurcation 
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while the threshold intensity, determined by dX_ /dip — 0, is 

/t,, = /o + AiA2 

which is minimum for Ai/ki — A2/K2 and equal to K1K2 + A1A2, like for the conventional DRO. 
The signal intensity I2 is related to Ii via the relation 

iai/a2f = A2/A1 



(18) 



(19) 



deduced from Eqs.(14a-b), so that the two detunings must have the same sign. 

The solutions and X_ are represented by the lines a and a', respectively, in Fig. 2 that displays a subcritical 
bifurcation: Indeed, only the stationary solution X_^_ is stable ( solid line a ) and extends from the saddle-node intensity 
Iq to the Hopf bifurcation intensity Ih, above which the solution is periodic (dashed portion) . The stationary solution 
X is marginally unstable (dashed line a'), whatever the input intensity might be [See Appendix]. 

Let us notice that the conventional DRO detuning condition, A= Ai/ki = A2/K2 does no longer necessarily hold 
and hence the SPL-DRO can oscillate with a wider detuning range. This result is similar to the result of a 2:1 
degenerate SPL-DRO, in the case of a linear coupling. However, this latter case does not display any subcriticality: 
The linear coupling gives rise instead to two self phase-locked states corresponding to two distinct thresholds, for a 
given pumping rate |£1[ . 

The phases fj, and rj fulfill the relations 



with 



Qfp cos /i - 



sm /i 

sin 7] 
\x\ ai COST] 



1 

2a 



p 



1 




H2 





K2 




-VA1A2 



(20a) 

(20b) 
(20c) 



Eqs.(20) completely determine the absolute phases of the subharmonic waves, which are hence self-locked. These 
relations on the idler and signal phases, while similar to those in ref will be shown to present substantial 
differences. As a consequence of the self-locking, highly phase-coherent subharmonic outputs are expected from the 
SPL-DRO. 



In the case Ai_2 = 0, the bifurcation is supercritical with 

X = Ip ~ Ith 
hh = K1K2 

The real and imaginary parts of Eqs. (14a-b) give rise to 



(21a) 
(21b) 



sm fj. 



sm?7 = 



2ar, 



ai a2 

Ki h K2 

a2 ai 



1 



2^/p - K1K2 



ai a2 

Kl K2 

a2 ai 



\p cos /i + \/Ip — K1K2 cos 77 = 



(22a) 

(22b) 
(22c) 



Eq.(22c) implies that cos/i and cosry are either of opposite signs or simultaneously nil. The first case corresponds 



to 



— = — = 1 with sm 7] = v and sm ^ = 

a2 K2 Up 



(23) 



The case cos = cos 77 = is compatible only with ki 7^ ^2 , leading to 
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sin 77 = ±1 and sin/i = — 1 



(24b) 



where N — Ip/ kiK2 is the pumping rate. Eq.(24a) displays a dependence of ^ on the pump amphtude, unhkc eq.(19) 
and consequently the solutions are not continuous when Ai 2 — > 0, in the general case Ki 7^ K2 • 

For A1.2 = 0, A2,i 7^ 0, Eq. (15) has no real solution. Actually the exact solutions will be shown to display a 
time-dependent regime. 

Finally let us note that the signal and idler intensities are inversely proportional to |x|^ as displayed by Eq.(16). 
This dependence is a consequence of the MLl approximation, Eqs.(9). (The limit \x\ — > is irrelevant since the 
meanfield equations (13) are misleading for the conventional DRO). 

B. Numerical results 

The time evolution of the field amplitudes is obtained by solving the propagation equations (5-6) with the boundary 
conditions (11) for a given input amplitude and small initial signal and idler amplitudes. The solutions are obtained 
by numerical integration of Eqs. (5-6) using of a fourth-order Runge Kutta algorithm, until convergence is achieved. 
In most of the calculations, the cavity loss coefficients are taken constant and equal, ki^2 = 0.005, and the SHG phase 
mismatch ^ = t^kLijl = 0, unless otherwise stated. The input amplitude, the detunings and the coupling constant 
are varied. 

• Stationary solutions 

As predicted by the linear stability analysis of the meanfield solutions, a single solution, corresponding to , occurs 
by numerical integration for an input amplitude above the threshold given by Eq.(18), where the trivial solution is 
unstable. In a general manner, the signal and idler intensities satisfy the relations predicted by the mean-field, cither 
for non-zero detunings in Eq. (19) or vanishing detunings in Eqs. (23-24). Nevertheless, for zero detunings, numerical 
intensities agree with the meanfield intensities, only for small coupling parameters (|x| 0.15). As already pointed 
out, there is no continuity for the intensities when Ai^2 — > in the case k\ 7^ K2- 

The exact stationary amplitude |x| ai and its stability domain depend on the magnitude of |x| as displayed in Fig. 
2 by line h for |x| = 0.05 and line c for |x| — 0.2. In this latter case, the exact solution is close to the meanfield one, 
the stationary solution extends almost from Jp in the subthreshold domain, reachable only by backward adiabatic 
decrementation of the pump intensity, to the Hopf bifurcation threshold intensity J™™ , slightly larger than the 
meanfield value. As |x| decreases, the domain of stability is shortened below threshold (/q"™ > /g ),but it is enlarged 
above threshold, because the Hopf bifurcation threshold is shifted towards larger input intensity, > ( For 

Ki 2 = Ai 2 = 0.005, the meanfield predicts Ih = 2/t;i, while the numerical values are l"^"^ ~ 2.6Ith for |x| = 0.2 and 
jnum > 4j^^ for Ixl = 0.05). 

Subcriticality is also evidenced in Fig. 3, when starting the numerical integration from A= and performing an 
adiabatic increase of A, for a fixed input pump intensity Ip — 4kiK2. The different curves in solid line correspond to 
different values of |x|- They are symmetric for negative detunings, so that the tuning curve of the SPL-DRO appears 
as a widened double-sided fringe, which reminds of resonantly phase-modulated Fabry Perot devices (the same kind 
of fringe is obtained from a resonator containing an electro-optic phase modulator driven by a RF oscillator whose 
frequency is equal to the resonator free spectral range ^^). In the case of the conventional DRO (|x| = 0, dashed 
curve), lasing begins at A= where the intensity is maximum as seen in Eq.(12c) and stops for A= \/3^ indeed, the 
bifurcation is supercritical so that lasing may occur only if the input pump intensity is larger than the threshold value 
K1K2 + A1A2. The case |x| — 0.01 (curve a), corresponding to nascent bistability, displays lasing, approximately in 
the same range of detuning as in the |xl = case, with the significant difference that there is self phase-locking (see 
below on Fig. 6). As |x| increases further, the bifurcation becomes subcritical, the saddle-node moves away from the 
threshold, approaching the meanfield location, independently of A. The intensity reaches the same maximum /max for 
any < |x| < 0.25, but at a detuning Aniax( |xl)i for which the pump is entirely depleted; then the intensity decreases 
to zero when increasing further the detuning. The maximum /max is equal to the maximum conventional DRO 
intensity, which occurs at Amax = 0. In case of the parameters used in Fig. 3, Eq. (12c) gives rise to V /max = 0.1, in 
agreement with the numerical result. The detuning Amax becomes approximately proportional to |x| for large enough 
couphng strength. For |x| > 0.25, the idler intensity becomes time-dependent before Amax is reached. 

In summary, the cascaded SHG nonlinearity induces a |xl ~ dependent phase that shifts the optimum detuning 

from to Amax • 



7 



The temporal response of the system subject to stepwise- Uke detuning jumps has been also studied in order to check 
the stability of the subthreshold states against external perturbations that tend to modulate the detuning parameter 
(via the cavity length for instance). Figure 4, associated with |x| = 0.2, shows that the system recovers its steady 
state operating point as long as the perturbation amplitude does not exceed ±0.12 Amax. Note the longer decay time, 
which is characteristic of a critical slowing-down phenomenon, on the positive detuning step side that brings the DRO 
very close to the saddle-node bifurcation. 

We can define a self-locking domain in the (Ai , A2) plane over which the DRO is self phase-locked with a well-defined 
phase relationship. Contour plots of the numerical solutions for the idler Ii, the signal I2, the total subharmonic 
intensity Ii + I2 and the phase difference rj are shown in Figs 5 for \x\ = 0.2 and for a given pumping rate A'' = 
Ip/KiK2 =4 . The subharmonic intensities are scaled to the conventional zero-detuning DRO values as deduced from 
Eqs.(12c), //'^'^ = 2k2,i- These contour plots are obtained by adiabatical following of the stationary solutions in 
order to reveal the 2D subthreshold domain. The granular small regions adjacent to the Ai = axis correspond to a 
time-dependent regime. The h and I2 intensity distributions in Figs 5a-5b display off-diagonal maxima that are larger 
than the conventional DRO signal and idler outputs ( I\^2/Ii2'^ = 1-6). Note also that the exact intensities I\ and 
I2 are not invariant with respect to the product A1A2, differently from the predicted meanfield intensities (16) but 
they satisfy the relation (19). The figure (5c) shows that the maximum total intensity, slightly larger (/i + 12 = 2.2 ) 
than the conventional DRO total intensity, Ii -\- 12 = 2, oc;curs in the subthreshold domain close to the saddle-node 
bifurcation (see diagonal line d in Fig. 3). These curves show that the relative phase control of the subharmonics can 
be achieved via the control of the relative output intensities, provided that an independent control of both cavity 
detunings is implemented (see section V). 

• Phase-locking 

The numerical phase 77 related to the difference between the signal and idler phases agrees extremely well with the 
mean- field prediction in Eqs. (20b), ( 22b) or (23), (24b) for any values of the detuning and cavity loss parameter, but 
the numerical phase fx related to the sum of the signal and idler phases is found to depart significantly from the mean- 
field value (20a) when Ai^2/ki,2 ^ 1- We have checked numerically the self phase-locking effect of the signal-idler 
output, by varing randomly the initial phase of the signal and idler noise. This general result is illustrated in Fig. 6, 
for Ai/ki = A2/K2 where the phase difference r] follows randomly the initial values (blank circles), when |x| = 0, as 
expected for a conventional DRO. However, yet for a vanishingly small value S = 0.001, the phase difference jj locks 
to or TT (mod kw) for ^ = and the sum phase locks to a constant value determined only by the pumping rate. 

• Time-dependent solutions 

For Ai,2 7^ and Ai/ki = A2/k,2, the periodic solution remains generally stable on a large range of the input 

intensity. For instance, for |x| = 0.2 and A 1 (case of line c in Fig. 2), the periodic solution is stable until the 
pump parameter reaches OLp/octh ~ 12 , well above any experimentally achievable input. The period of the amplitude 
oscillation in the vicinity of 7^""* is found to be T ~ 300 r. But the case |Ai/ki — A2/K2I 7^ leads to a more 
complex dynamics, which depends on the nonlinear coupling \x\ and departs from the meanfield predictions . 

In the zero-detuning case, no time-dependent solutions arc found for realistic pumping rate {N < 25) for equal 
cavity loss ki = K2- Differently, the case ki ^ K2 displays time-dependent regimes, for \x\ ^ 0.2. In the case 
K\l K2 = 2.5, Ixl = 0.2 shown in Fig.7, a periodic regime arises at /j„ = 2 liy^ and a weak chaotic pulsating behaviour 
is found at /j„ = 4 Jt/j, with a pulsating signal intensity. Vanishingly small detunings (Ai^2) however give rise to 
steady state solutions satisfying relations (16). 

The numerical solutions associated with Ai^2 = 0, A2,i 0, are not stationary, as predicted by the meanfield 
model; they are periodic with time for any input above threshold. 

• Influence of SHG phase-mismatch 

We have studied the influence of a moderate SHG phase mismatch (^ ^ 0) on the dynamical behaviour of the 
system. This phase mismatch is usually controlled by the temperature or the angular orientation of the SHG crystal. 
It provides a control of the strength of the nonlinear coupling parameter |x| , via the relation (10), and offsets the 
value to which the phase difference t] locks (Fig. 6). For instance, a small phase ^ < 1 does not significantly change the 
idler amplitudes of Fig. 3 but induces a slight imbalance of the signal and idler intensity ratio compared to the ratio 
(19) and slightly reduces the self-locking range. The relation /i = (mod tt) is still valid and 77 locks to ^. In some 
undesirable operating conditions, for instance, when the system approaches a Hopf bifurcation for a given pumping 
rate and detuning set (see Fig.2 for S = 0.2, ^ = 0), a small amount of phase mismatch produces a shift of the periodic 
oscillation threshold towards higher pump rate. Differently, for A= and the same other parameters as for curve (c) 
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of Fig. 2, a small phase mismatch = 0.1) leads to a slow periodic regime (T ~ 4 x lO"*"^ r at ain/ath = 1-2) which 
is suppressed when vanishingly small detunings (always present in practical devices) are introduced. 

From the simulations carried out, we conclude that a moderate SHG phase- mismatch does not change significantly 
the main bifurcation dynamics studied up to this point, except for the Ai_2 = case. 

Finally let us point out that the SPL-DRO solutions expanded up to the second-order in Eq. (9) (ML2 approxima- 
tion) agree very well with the exact solutions in almost the whole range of the pump intensity and detunings below 
the Hopf bifurcation threshold, except for the singular case Ai^2 — corresponding to the transition from subcriti- 
cality to supercriticality, where it fails to converge. Only the propagation model can solve this case, which requires a 
double-precision computation due to the slow convergence associated to a critical slowing-down phenomenon. 



IV. SPL-TRO CASE 



In the case of SPL-TRO, the meanfield equations (13b)-(13c) are still valid, except that /p = denotes now the 
circulating intracavity pump intensity. The intracavity pump amplitude obeys the boundary condition (11a), leading 
to 

[1 - Tp exp(iAp)]Ap - iTp exp(iAp)yliyl2 = A„ (25) 

where the pump detuning may be arbitrarily large (Ap < 27r) . The relation between the input pump parameter 
lin — \Ain\^ — Q;f„, and the cavity pump intensity is obtained, for an arbitrary phase of the input field 

Im = (1 + 7'p - 2rpCos Ap)Q!p + r^a\a\ + 2rpaia2ap[sin(Ap - ^i) + rp sin/i] (26) 

where ai,a2 and (i can be deduced from Eqs.(16) and (19-20). The threshold for oscillation is easily found from 
Eq. (26), on the basis of simple considerations. Below the threshold, where ai_2 = 0, Ip grows linearly as a function 
of lin- When it reaches the value given by Eq.(21), oscillation starts, which leads to the following threshold input 
intensity 

jSPL-TBO ^ (1 + ^2 _ 2rpCOS Ap)(/o + A1A2) (27) 

In a practical TROs, although the pump finesse is lower than the finesse at the subharmonic waves (typically an 
order of magnitude), the condition Kp ^ 1 is still satisfied . If we consider a small enough pump detuning and the 
specific case A= Ai/ki — A2/K2 , Eq.(27) reduces to the conventional TRO threshold js^ 

/™ = \A,n\l = 4^i^2{l+ AK1+ Ap) (28) 

where Ap= Ap/np. Fig. 8 shows the normalized intra-cavity pump and idler bifurcation diagrams as a function 

of \/lin = \J lin/ljf^^ when condition (12a) holds and for Ap = 0. Only the details of the \x\oii solutions in the 
vicinity of the threshold are plotted in the top frame of Fig. 8 that displays the subcriticality. The numerical and 
meanfield solutions of the idler amplitude as given by Eq. (16) are confounded over a much larger input intensity 
range. The thick solid and dashed lines for ap correspond to the meanfield solution (26), while the thin solid line is 
the exact numerical pump solution. Notice that the numerical solution agrees very well with the meanfield solution. 
As the input parameter is increased from threshold, Ip decreases from the clamped value of the conventional TRO to 

a minimum /q (Eq.17) for an input pump intensity /j„ given by 

2 _2 

/(1+A) (29) 

for Ap = 0. The meanfield model hence predicts that, unlike for conventional TRO, the intracavity pump intensity is 

not clamped. For A^ 0, /j„ 1, e.g. there is a transition to supercriticality as in the SPL-DRO case. 

In pump-detuned conventional TROs, subcriticality occurs only when the condition ApA> 1 holds [Q. The 
experimental observation of this subcriticality requires a high pumping level because of the pump detuning dependence 
in Eq.(28) ||3^ . It is interesting to study how the subcriticality originating from the nonlinear OPO/SHG cascading 
would affect the intrinsic detuned TRO bistability curve. In Fig. 9 we have plotted the numerically computed idler 
intensity versus the normalized input intensity for \x\ = 0,0.2,0.6 and ApA= 2 . The subthreshold domain extends 



-mm 

lin. 



l + («iAp)(A/|xl)' 
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as Ixl increases from zero, it has also increased, when compared to the Ap= 0, |x| = 0.2 case, shown in Fig. 8. 
Furthermore, as |x| increases, the numerical threshold is smaller than the meanfield threshold given by Eq.(27). 

We have also verified that the numerical amplitudes and phases agree with the meanfield predictions for 7i„< 10 
and for detuning values such that Aj/k^ ;^ 1. Far above threshold and for larger detunings, only the phase difference 
relations remain exact, like in the SPL-DRO case. 

Furthermore, the detuning range A for self-locking is less dependent on the coupling strength than in the SPL-DRO 
case. 

The high-pump-fincssc SPL-TRO device does not exhibit the time-dependent periodical solutions observed above 
threshold with the SPL-DRO as long as |x| < 0.52 (which corresponds roughly to equal OPO and SHG crystal length 

for a PPLN) and for realistic maximum pumping rate Im— hn/ Ith^~"'^^'~' < 2500, corresponding to a Watt-level 
pump power and typical TRO pump thresholds in the mW range. For larger coupling |x| = 0.6, periodic oscillations 
occur, with a period varying in a very complex way when the input pump is increased. Nevertheless no chaotic regime 
is observed in the range of considered input pump. 

In many practical experiments, one cannot avoid a weak pump resonance in DROs due to the multiwavelength 
coatings of the mirrors. It is thus interesting to investigate how a moderate pump resonance would affect the (Ai, A2) 
self-locking domain of Figs. 5. The input external pump intensity is kept equal to /f^* — 4kiK2 to provide a comparison 

with Figs. 5 ( Ixl = 0.2). Such a pumping level corresponds to an internal TRO pumping rate /j„ Z^™*^ — 8/Kp = 160. 
In Fig. 10, the pump reflectivity is rp = 0.8 with Ap = 0. Apart from the slight enlargement of the self-locking 
domains, as compared to Fig. 5, which is mainly attributed to the high pumping rate, the intensities are smaller than 
for the DRO, also the intensity distributions are strongly modified, with off-diagonal maxima. 

In summary, the only improvement due to a moderate pump resonance is the stabilization effect respective to the 
onset of temporal dynamics. The extended self-phase locking range is paid back with lower output intensities. 

V. PRACTICAL IMPLEMENTATION OF SPL-OPOS 

In order to avoid spurious cavity loss, the use of a dual-grating quasi-phasc matched periodically poled crystal 
is particularly well suited to the implementation of SPL-OPOs. However, the grating period should be accurately 
designed so as to phase-matched simultaneously both interactions for the same chip temperature, even though the 
theoretical analysis predicts a minor influence of SHG phase mismatch. 

From an experimental point of view, it is desirable to implement 3:1 SPL-OPOs using the least constraining OPO 
configuration. Diagnosis methods to check the high phase coherence between the subharmonics under SPL have to be 
implemented. In the frequency domain, this can be achieved by monitoring the beatnote between the signal wave and 
the externally frequency-doubled idler wave (or between the summed subharmonics and the pump). Because under 
SPL operation these two waves are frequency degenerate, the output signal wave must be preliminary frequency- 
shifted by a suitable RF frequency unp using, e.g., an acousto-optic modulator. When the OPO operates within the 
locking range, the beatnote frequency should be fixed to unp and its power spectral density should approach a Dirac 
function. Another equivalent method to check the phase coherence is to perform an interferometric fringe pattern 
measurement by overlapping the two beams on a slow detector ||l7|| . In the following, other indirect methods, based 
on the theoretical analysis will be outlined. 

As a starting point of the analysis, we have assumed perfect 3:2:1 frequency ratios for the pump, signal and idler 
waves. In practical devices, such a situation will be unlikely met at once. The major difficulty will come from the 
fine (continuous) tuning of the signal and idler frequencies close enough to the 2:1 degeneracy, in order to fall within 
the capture range of the self-phase locking. For a fixed pump frequency, such a fine tuning is usually performed 
via the temperature or angle tuning of the phase-matching. Singly resonant devices (PRSROs) offer an easy and 
relaxed mode- hop free frequency tuning because only one subhamonic wave is resonant, especially in a dual- arm 
cavity configuration to control independently the pump and signal resonances p^ ] . However from our analysis of the 
SPL-PRSRO with a resonant signal wave, which also displays subcriticality, extremely small self-locking ranges are 
predicted due to the weak coupling between the injecting frequency-doubled idler and the resonant signal. While in 
a conventional PRSRO the signal field is constrained to oscillate with a nil cavity detuning, the nonlinear coupling 
is found to allow oscillation over a small detuning range not exceeding the cavity linewidth. The increase of the 
locking range versus the coupling parameter (e.g. the SHG crystal length) is only moderate, even with |x| — 1 (which 
would correspond to a SHG crystal ~ 1.5 times longer than the OPO crystal supposing that both OPO and SHG 
interactions have the same nonlinearity magnitude) . Due to this limited capture range, experimentally confirmed in 
ref 0, 3:1 SPL-PRSROs would probably require an additional electronic servo on the cavity length to operate as a 
stable divider device. A convenient criterion for the assessment of SPL in PRSROs would be the slight enhancement 
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of the signal intensity when the nonhnear couphng is switched on, compared to the sHght decrease predicted for the 
SPL-DRO (Fig. 3 ). A possible way to realize such a switching is to have an OPO-only grating section patterned 
beside the OPO/SHG dual-grating section on the periodically-poled wafer. 

From our theoretical study, the widest SPL range (a few cavity linewidths) is obtained with the doubly resonant 
configuration (SPL-DRO) due to the strong self-injection regime, with eventually a weak pump resonance (SPL-TRO) 
to stabilize the device. Even though the amount of nonlinear coupling required can be extremely small ( |x| = 0.01 ) 
- and such low level doubled idler can even be spontaneously generated via non-phase matched or higher-order quasi- 
phase matching in PP single-grating nonlinear OPO crystals - a coupling strength corresponding to |x| = 0.1 — 0.2 will 
ensure a robust self-phase locking of the subharmonic wave. Preliminary single mode-pair operation of DRO /TROs 
usually require an intensity sidelock servo to control the stability of the cavity length. This sidelock servo compares 
the output signal (or idler) intensity to a stable electronic voltage reference which sets the operating (usually non-zero) 
signal and idler cavity detuning values. When a linear cavity is used, the signal and idler (plus eventually the pump) 
detunings cannot be independently controlled via the cavity length. It is then probable that the oscillating mode pair 
will have equal normalized detunings (condition (12a)) that satisfy the minimum threshold. Under sidelock servo the 
transition from conventional to SPL states is accompanied necessarily by a detuning transition (see Fig. 3). It is then 
important to set the sidelock reference voltage as close as to the maximum fringe intensity in order that the new 
detuning value does not exceed the allowed subthreshold range. In the case of a well resolved DRO mode pair cluster, 
the observation of these subthreshold states, and the associated broadened mode pair fringe, should be made possible 
via adiabatic cavity length tuning. Such an observation would be an indirect diagnosis of SPL. But it is necessary to 
have an independent control of the signal and idler detunings to explore the full allowed range of (Ai, A2) detunings 
depicted in Figs. 5 and 10. A dual-cavity DRO/TRO design would then be appropriate. The control of these detunings 
allows the control of the relative phase between the pump, signal and idler. The output subharmonic intensity ratio 
/1//2 (see Eq.l9) can be used as an error signal for the relative phase control. 



VI. CONCLUSIONS AND OUTLOOK 



We have theoretically demonstrated that resonant x*^^-' : x*^^' nonlinear cascaded OPO/SHG processes induce a self 
injection-locking between the subharmonic waves of an OPO leading to the self-phase locking of the three interacting 
waves, unlike in a conventional OPO for which the absolute phases of the signal and idler are undetermined. The 
theoretical treatment encompasses the meanfield model as well as a full propagation model. The doubly and triply 
resonant oscillator configurations lead to the widest self phase-locking ranges. The main conclusions are: a) the 
minimum threshold of oscillation of these devices are identical to the conventional devices; b) the nonlinear cascading 
leads to a subcritical behaviour, even in the case of a DRO or a PR-SRO, and can lead to the occurrence of self-pulsing 
instabilitities. This subcriticality is different from the standard subcriticality reported in pump detuned conventional 
TROs. The nonlinear coupling removes the detuning constraints of conventional systems, allowing for a potentially 
accurate control of the relative phase between the subharmonic waves. The range of allowable detunings over which 
the field phases are locked depends on the magnitude of the nonlinear coupling and on the self-injection regime. 
While this range is smaller than the cavity linewidth for PRSROs (weak injection regime), it spans over several 
cavity linewidths under strong injection-locking regime obtained in signal/idler resonant devices (DRO/TROs). The 
SPL-DRO/TRO give rise to a richer dynamics than the singly resonant PRSRO for which no occurence of a Hopf 
bifurcation is found in the whole range of the system parameters. It is thus interesting to extend the theoretical 
model by including diffraction effects in order to investigate the possible occurrence of new spatio-temporal dynamics, 
p3[ , |37|. These self-phase locked OPOs will be useful tools for applications requiring a high degree of optical phase 
coherence between optical harmonic waves, such as precision optical measurements in the mid-IR or Fourier synthesis 
of ultra-short optical pulses. 

The model developed can be easily extended to the study of divide-by-4 SPL-OPOs based on the cascading 
OPO/OPO processes Aoj ^ 2uj ^ u which has the potential to generate up to 8 phase- locked harmonic waves 
by additional up-conversion processes. Such a strong nonlinear ly coupled system can be viewed as a secondary degen- 
erate OPO (DOPO) embedded in a primary DOPO. Our future work will be directed to the theoretical investigation 
of the stability of 4:1 OPO dividers. A classical signature of the system, derived from the meanfield analysis, is the 
clamping of the secondary pump 2lu to the threshold power for the fundamental oscillation. 

The present study of 3:1 OPO dividers, which makes the simplest assumption of exact 3:2:1 frequency ratios, arouses 
another interrogation. An interesting situation not considered regards the behaviour of the nearly 3:1 OPO/SHG 
system when the frequency ratios slightly departs from the perfect 3:2:1 division ratio by a radio-frequency quantity 
S « uj, e.g. when Bui 2iLj ~ S,lj + 6. The frequency difference beween the signal wave and the doubled idler is 
then \2uji — uJs\ = 35. When the doubled idler frequency does not match one of the cavity eigenmode frequencies, the 
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operation of the OPO would be merely that of a conventional DRO. But if the OPO is tuned such that 3d = FSRg ~ 
c/A, FSRs being the signal free spectral range of the cavity, then the doubled idler will be enhanced to a point where 
it may lead to the creation of a new mode pair with frequencies {2uj + 2S,uj ~ 26) and so on. The parametric gain 
bandwidth of OPOs extending usually over several THz or several tens of THz (in case of a wavelength non-critical 
phase matching), a multitude of self- phase locked mode pairs equally spaced by 3d may potentially oscillate, provided 
that the pumping level is sufficiently high. Such a complex system opens the prospect of building a mode-locked, 
dual-band OPO frequency comb generator using cascaded second-order nonlinearities as the passive mode-locking 
mechanism. In the time domain, the output of such a system would consist of a train of short optical pulses with 
a repetion rate set by the FSRs intermode spacing, provided that the relative phase between adjacent mode pairs 
is preserved and group velocity dispersion is compensated. We note that a similar cw-DRO running near frequency 
degeneracy with thousands of mode pairs actively locked by an intracavity electro-optic phase modulator has been 
recently reported, with a striking passive output stability feature compared to a conventional quasi-degenerate single- 
mode pair DRO [|§. 
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n° ERBIC 15CT 98 0814 ). 



VII. 



APPENDIX I: LINEAR STABILITY ANALYSIS OF THE SPL-DRO 



With Aj{t) = Aj{t,0), Eqs. (9-11) give rise to the mapping equations 

A2{t + r) = r^e'^^lA.it) + iAp{t)A;it) + ix*Alit)]] 
A,{t + r) = r,e'^-[A,it) + iAp{t)Alit) + ixA^mUt)] 

Ap{t) — Ain 



(Al) 
(A2) 
(A3) 



the stationary solutions of which are deduced for large ri_2 and small A1.2. (See Eqs.(19),(22)-(23)). 

The linear stability analysis consists in assuming small deviations dAi^2{t) from the stationary solutions Ai_2 , 



^l,2(i) =^l,2 + Mi,2(0 



(A4a) 



X 



(A4b) 



where A may be complex. The stationary solutions Ai_2 are stable only if the real part of any A is negative. At 
the instability threshold, therefore the system may undergo a Hopf bifurcation ( Ao = ±i/?), so that the intensitiy 
oscillates with time at angular frequency /3. 

Eqs. (A1-A4) lead to a linearized system of four equations, the determinant of which satisfies 



D = 



A + ai 
^XIA^_ 
-2ix*2Ai 
iA*„ o 



-«Xi^2 -ixiAi -iXiAins 

A + al iAl^ ixlA, 

~iAijn,2 A + 0:2 

2%^! A + a*2 



= 



(A5) 



with the notations 



A: 



(A6) 



and 



ai,2 = Kl,2 - «Ai_2 

Xi = e'^'xri,X2 



Ai' 



,-'iAi 



im,l,2 = '^1,26 



X^2 



(A7a) 
(A7b) 
(A7c) 
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Then, the eigenvalues A are solutions of the quartic characteristic equation 



A^ + $3A^ + $2A^ + $iA + $0 = (A8) 

where all the coefficients $j are real 

$3 = -2(ki + K2) (A9a) 

2 2 I 1 2 I 1 2 

*2 = |ai| +\a2\ + 4kiK2 - 2rir2 |Ai„| cos(Ai - A2) - |xr2A2 1 

+4rir2 |x^i 1^ cos(Ai + A2) (A9b) 
$1 = 2rir2 |A,„f [(Ai + A2) sin(Ai + A2) - (K1+K2) cos(Ai - A2)] + 2(^2 bil' + ki |a2|') 

-2fi;2 1x^2^2^ +4|xrir23if + K2) cos(Ai + A2) + (Ai+ A2) sin(Ai + A2)] 

~inr2 \x\^ cos A2(Ai„AlAl - c.c.) + Srira \xf sin A2{AnAlAl + c.c.) (A9c) 
$0 = |ai|^ |a2|^ - 2rir2 \An\'^ [(kiK2 + A1A2) cos(Ai - A2) + (A1K2 - A2K1) sin(Ai - A2)] 

+4rir2 1x^1 r [(«i'^'2 - Ai A2) cos(Ai + A2) + (A1K2 + A2K1) sin(Ai + A2)] 

- \xria2A2\^ + (4|xr2Aif -r^ |A,„f) (jxri^if -r? pi„|^) (A9d) 

The roots for A and consequently the eigenvalues A have been calculated for various detunings Ai^2 as a function of 
the input amplitude Ain for the trivial solutions Ai^2 = and the non trivial solutions {Ai^2)± deduced from Eqs.(16) 
and (20). 

The trivial solution is verified to be stable for any input pump intensity below /t/((Eq.l8). 

On the other hand the sohition X is found to be unstable with respect to any constant perturbation for any 
detuning and input pump intensity: Indeed, in this case, Eq.( A8) has a root A = 0,i.e, A = 0, because the constant 
coefficient $0 is identically nil, when replacing X_ by its expression (16). 

Differently, the upper branch displays two sets of complex conjugate eigenvalues, (cri,2 ± 1^1,2) with negative real 
parts for a range of the input pump intensity, lying from Iq to Ih where the system undergoes a Hopf bifurcation 
((71 = 0). For Ip > Ih, the signal and idler intensities oscillate with time at angular frequency Therefore the 
upper branch is stable for Io<Ip < Ih- 

The Hopf bifurcation threshold intensity Ih may be either smaller or larger than the threshold intensity Ith, 
depending on detunings Ai and A2. Let us introduce the parameter fxn that measures the departure between the 
threshold for lasing 1^^ and the Hopf bifurcation threshold as 

M/f = {Ith - lH)/Ith (AlO) 

The variation of ij,h is presented in Fig. 11 as a function Ai = A either with Ai = A2 in curve (a) or Ai = iA2 
in curve (b). 

For Ai = A2, the Hofp bifurcation (cti = 0,(T2 < 0) occurs above threshold only for A < k. Therefore lasing 

stationary solutions can be reached from rest only for detunings smaller than the cavity loss coefficient. Otherwise, 
lasing at larger detunings can be reached, when increasing adiabatically the detunings from values smaller than k. In 
the other case, Ai — ^^2, i^h is positive for any Ai < 5k. For higher detuning, a2 reaches zero at an input intensity 
much smaller than the value at which cti = crosses zero. This causes an abrupt decrease of Ih, so that fxn suddenly 
becomes negative for Ai ^ 5k, as shown in the curve (b) of Fig. 11. 

Finally, the angular frequency (3h at the Hopf bifurcation threshold, not reported here, is found to vary propor- 
tionally to the detuning. 
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Figure captions 



Fig.l: Schematic ring cavity model of SPL-OPOs. All intracavity losses are lumped into the output mirror trans- 
missivities tj. 

Fig. 2: Comparison of the mcanficld and numerical stationary solutions of SPL-DRO, with ki^2 = ^1,2 = 0.005, 
as a function of the scaled pump input amplitude {ath = Vhh)- (a): meanfield upper branch (The dotted portion 
corresponds to unstable solutions), (a'): meanfield lower branch . (b): numerical solution for |x| =0.05. (c): numerical 
solutions for |x| = 0.2. 

Fig.3:Numerical SPL-DRO self locking ranges as a function of the detuning Ai/ki(= A2/K2) , for Im = 4kiK2, with 
Ixl = 0.01 in (a), 0.05 in (b), 0.1 in (c), 0.2 in (d) and 0.25 in (e). The dashed line is for |x| = 0. The parameters are 
2 = 0.005. 

Fig. 4: Time response of the SPL-DRO stationary state operating at = 0.2, Aopt= 4|, imder step-wise linear 

detuning jumps Kopt with modulation index /? = 0.109 (curve a). Curve (b) gives the idler amplitude response 

and curve (c) the phase difference response. The re-capture range of the detuning perturbation is /3max = 0.12. 
Fig. 5: Contour plots of the stationary SPL-DRO signal (a), idler (b), and total (c) intensities and phase difference r] 
(d) in the 2D detunig plane, obtained from the numerical computation with ki,2 = 0.005, |xl = "S* = 0.2, /i„ = Akik^ 
. The intensities are scaled to Ji,2 = 2k2,i- The pale blue domain are trivial solutions. 

Fig.6: Numerically computed distribution of stationary signal/idler phase difference as defined in eq.(14d) versus 
the initial random phase for Ai^2 = 1^1,2 = 0.005, am/ ath = 2 and S = 0.001, for a SHG phase mismatch ^ = n (i.e 
X = 0) in blank circles; ^ = in solid black circles and ^ = 37r/2 in black triangles. Note that each phase data point 
is associated to the same stationary signal-idler intensities and the same phase sum /z. 

Fig. 7: Time-dependent numerical solutions for ki = 0.005, K2 = 0.002, Ai_2 = 0, x = 0.2, 7j„ = 1.96«;ik;2 (bottom 

frame) and = 4kiK2 (top frame). 

Figs. 8: Meanfield and numerical SPL-TRO stationary solutions for Ap = 0, = 0.9 and the same other parameters 
as for Fig. 2. The top frame shows details of the hysteresis loop of the idler amplitude in the vicinity of the input 
pump threshold. The thin lines correspond to the numerical solutions, interrupted by the vertical thin dashed line. 
The thick lines correspond to the meanfield solutions. The thick dashed line in the bottom frame is the meanfield 
unstable branch of ap/ath {ath = \/i^\i^2 + A1A2), and the thinner horizontal dashed line shows the conventional 
TRO pump clamping. 

Fig.9: Subcriticality in the strongly pump-detuned conventional TRO {S = 0) and SPL-TRO [S = 0.2, 0.6) for 
Tp = 0.9, A= 1, Ap= 2 and same other parameters as in Fig. 2. The thresholds located at the vertical arrows, are 
lower than the one given by eq.(27). 

Fig. 10: Contour plots of the numerical stationary SPL-TRO signal (a), idler (b), and total (c) intensities and phase 
difference rj (d) in the 2D detuning plane with Vp = 0.8, Ap = 0, Ki,2 = 0.005, x = 'S' = 0.2, Ji„ = 4k;ik;2 . The 
intensities are scaled to /i,2 = 2k;2,i. In the pale blue domain there is no lasing. For = 0.9, the self locking domain 
extends over the whole 2D frame. 

Fig. 11: Plots of the reduced Hopf bifurcation thresholds /i/f = {Ith — lH)/Ith versus the normalised idler detiming 
Ai/ki, for A2 = Ai (solid line (a)); A2 = 2Ai (solid line (b)). The cavity loss parameters used are Ki^2 = h = 0.005 
and S = 0.2. 
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